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Abstract 

A new version of the semiempirical formula based on fission approach of alpha decay is derived, by 
using the optimum values of the fitting parameters determined for even-even nuclei, combined with 
hindrance factors for even-odd, odd-even, and odd-odd nuclides. The deviations from experimental 
data for two regions of nuclear chart (493 alpha emitters with Z = 52 — 118 and 142 transuranium 
nuclei including superheavies {Z = 92 — 118), respectively) are compared with those obtained by 
using the universal curve and the Viola-Seaborg semiempirical relationship. 
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The existence of a fission barrier produced by shell effects in the region of superheavy 
nuclei continues to be proved by the successful synthesis of the heaviest nuclides yl 
which stimulate a corresponding theoretical development (see 
therein). In 1969 it was shown that fission barriers as high as 8 to 12 MeV develop due 
to the shell corrections around double shell closures which appeared in the standard shell 
model at Z = 114, = 184. Nowadays some other magic numbers have been proposed, 
e.g. Z = 126. The majority of these proton-rich nuclei are mainly decaying by a-particle 
emission, motivat ing an important effort to derive different methods allowing to estimate 
the half-lives 



tivat ing an important effort to derive different methods allo wing to esti 
|l^[ll[| and to perform systematic studies in this region e.g. [iJQ, 13]. 



Three such methods based on fission theory extended to a very large asymmetry 
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analytical superasymmetric fission (ASAF) mode^ 
and the semiempirical relationship (SemFIS) [ " 
The kinetic energy of the emitted particle, = QAi/A, is given at the input. It was 
shown that the description of data in the vicinity of the magic proton and neutron numbers, 
where the errors of the other relationships are large, was improved by introducing the SemFIS 
formula, in which the action integral is multiplied by a second order polynomial, x = xi^i y)i 
in the variables x and y expressing the "distance" from the nearest magic plus one neighbours 
x = {N- Ni)/{Ni+i - Ni) and y = {Z - Zi)/{Zi+i - Z,), where Ni > N > iV,+i and > 
Z > Zi^i. In this way not only the Z dependence (present in the majority of semiempirical 
relationships) but also the strong influence of the neutron shell effects are taken into account. 
A computer program allows to change automatically the fit parameters, every time 
a better set of experimental data is available. This polynomial approximation needs six 
parameters in each group of even-even, even-odd, odd-even and odd-odd alpha emitters. 
The purpose of the present work is to derive an alternative with less number of parameters 
and to compare the results with those obtained by using the universal curve. The same 
six quantities determined for even-even nuclei should be used in the other groups, only 
completed in the expression of log T by a different hindrance factor. 

We study a spontaneous process of a-decay in which a parent nucleus A, Z in its ground 
state is split into two fragments Ai = A — 4, Zi = Z — 2 and A2 = 4, Z2 = 2 

AZ ^'Zi+ ^''Z2 (1) 

in a way that conserves the hadron numbers A = A1 + A2, Z = Z1 + Z2. The heavy fragment 
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Ai, Z\ is called daughter and the light one emitted a-particle. The partial decay half-life T 
of the parent nucleus is related to the disintegration constant A of the exponential decay law 
in time by the relationship 

i^^A (2) 

The quantum mechanical tunnelling process leads to a relatively simple expression for A as 
the product of the three model dependent quantities 

A = uSP (3) 

where u is the frequency of assaults, S is the preformation probability of the a-particle at 
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the nuclear surface, and P is the penetrability of the external part of the barrier 

Very frequently the penetrability is calculated by using the one-dimensional Wentzel- 
Kramers-Brillouin (WKB) approximation 

Rb 



2 

I 

Ra 



P = exp(-J^); K = -j ^j2B{R)E{R)dR (4) 



where B is the nuclear inertia, approximated by the reduced mass = mAiA2/ A, m is the 
nucleon mass, E is the potential energy from which the Q- value has been subtracted out, 
Ra and Rb are the classical turning points, and K is the action integral. 

The equations |21ll] are used as a starting point by the three methods mentioned above 
(ASAF, UNIV, and SemFIS) based on fission theory of a-decay. The zero point vibration 
energy is Ei, = hv/2 in which h is the Plank constant. The outer potential barrier is of 
Coulomb nature, E{R) = e^ZiZ2/ R — Q, with the electron charge e = 1.43998 MeV-fm. 
The other numerical constants are given by (l/2)/iln2 = 1.4333-10^^^ MeV-s and 2\^2m/h 
= 0.43921 MeV-i/2fjn-i/2_ 

It is instructive to outline the derivation of the UNIV formula. For a-decay of even-even 
nuclei the eqs. |21ll] allow us to obtain a simple relationship of the decimal logarithm of the 
half-hfe 

logT= - log P + Cee (5) 

where the additive constant Cee = log^Q, — logz/ + log(ln2) = —20.325 if we are making two 
approximations: Sa = 0.0180302 and u = lO^^ *^^ s~^. For even-odd, odd-even and odd-odd 
nuclei we replace Cee by Ceo = Cee + h^o, Coe = Cee + hoe, and Coo = Cee + hoo, respectively, 
where heo, hoe, hoo are the mean values of the hindrance factors in these groups of nuclides. 



In a doubly logarithmic scale the equation El represents a straight line with a slope equal to 
unity. The penetrability of an external Coulomb barrier, having as the first turning point 
the separation distance at the touching configuration Ra = Rt = Ri + R2 and the second 
one defined by Z1Z2/ Rt = Q, niay be found analytically as 



logP = 0.22873(/iA^i^2i?b)^/^ arccos - Jr(l - r) 



(6) 



where r = Rt/Rb, Rt = 1.2249{AI^^ + A]^^), and Rb = IA3998Z1Z2/Q. We use the liquid 
drop model radius constant tq = 1.2249 fm. 
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A great number of alpha emitters are known 
figure n where we plotted the experimental half-lives of 162 even-even, 122 even-odd, 115 
odd-even, and of 94 odd-odd nuclides versus the neutron number of the daughter nucleus. 
The very strong neutron shell effect is clearly seen at the magic daughter number A^^^ = 126, 
where the shortest half-lives have been measured. 

The corresponding universal curves for a-decay of the same nuclides as in figure Q in 
four groups of even-even, even-odd, odd-even and odd-odd nuclei are displayed in figure |21 
The following hindrance factors have been used: hoe = 0.445, /leo = 0.294, hoo = 0.842, for 
even-odd, odd-even, and odd-odd nuclei, respectively. Qualitatively one can see that the 
data for even-even nuclei are well described by the universal curve. We can evaluate the 
standard root-mean-square (rms) deviation of logT values: 

( n ^ 1/2 

^=\jyog{T,/Te.,)]V{n-l)j (7) 

leading to cr = 0.352,0.612,0.546,0.841 for the 162 even-even, 122 even-odd, 115 odd-even, 
and of 94 odd-odd nuclides, respectively. 

Despite the fact that a main part of the strong shell effect present in the Q-value was 
accounted for, there are still some systematic errors in the vicinity of the magic number of 
ueutrous as may be seen in the plot of logT - logT^, against the neutron number of the 
daughter (see figure 14 of the ref. [20|). 

The same kind of underestimation around neutron magic numbers is also present in the 
semiempirical formulae. A typical example is the Viola-Seaborg [2^ semiempirical relation- 
ship 

log T={aZ- h)Q~^/^ - [cZ + d) (8) 
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We shall use for Z < 82, < 126 the parameter values a = 2.42151: 6 = 62.3848; c = 
0.59015, d = 4.2109, and for Z > 82, N > 126 the recently published [ll| ones a = 1.3892; 
b = —13.862; c = 0.1086, d = 41.458. It gives excellent agreement in the region of actinides 
but it underestimates the lifetimes of lighter nuclei. It is very frequently used, particularly 
in the region of superheavy nuclei. 

The behaviour around magic numbers can be improved by using the SemFIS formula in 
which we explicitly take into account the shell effect by multiplying the action integral, K, 
with the polynomial x = xi^yV) mentioned above: 

log T = 0.43429x(a:, y) ■ K - 20.446 + (9) 

where is a hindrance factor which takes different values if/g = —0.025 for even-even 
emitters, = 0.420 for e-o, H^^ = 0.280 for o-e, and = 0.810 for o-o ones. 

K = 2.52956Zi[Ai/(Ag)]^/2[arccosv^- yr(r^7)] ; 

r = 0.423g(1.5874 + A}/^)/Zi (10) 

The numerical coefficient x? close to unity, is a second-order polynomial 

X = Bi+ x{B2 + xBi) + y{B^ + yBo) + xyB^ (11) 

with coefficients obtained by fitting the data plotted in figure^ for 163 even-even emitters: 
Bi = 0.987389, B2 = 0.009520, B^ = 0.033711, B3 = 0.035249, Bq = -0.038941, B^ = 
0.019451. The reduced variables are defined 

x = {N- iV.)/(iV.+i -Ni); Ni<N< iV,+i (12) 

y=iZ- Z,)/(Z,+i - Z,) ; Z,<Z< (13) 

with Ni = ....,51,83,127,185,229, , Zi = ....,29,51,83,115, hence for the region of 

superheavy nuclei x = {Z - 83)/(127 - 83), y = {N - 127)/(185 - 127). 

The improvement can be seen by comparing the standard deviations of UNIV and SemFIS 
for the large set of data in table IH The figure El shows that there are nolonger large sys- 
tematical discrepancies around magic numbers, because shell effects were taken into account 
by SemFIS. The example are given for the best results obtained in the region of even-even 
(top) and odd-even (bottom) nuclei. 
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TABLE I: Standard deviations for semiempirical formula, and universal curves in the region of all 
alpha emitters (Z = 52 — 118). 



n parity 




CsemFIS 


162 e-e 


0.352 


0.233 


122 e-o 


0.612 


0.515 


115 o-e 


0.546 


0.497 


94 o-o 


0.841 


0.709 



In order to illustrate the possibility of using a different set of parameter values, we apply 
the semFIS formula for transuranium nuclei (the set of 47 even-even {Z = 92,118), 45 
even-odd {Z = 92,114), 25 odd-even {Z = 93,115), 25 odd-odd {Z = 93,113) nuclides). 
Now we obtain the following parameters: Bi = 0.985415, B2 = 0.102199, B4 = -0.832081, 
^3 = -0.024863, Be = -0.681221, B5 = 1.50572 and if£ = 0, = 0.63, H^^ = 0.51, 
HL = 1-26. 

The optimum hindrance factors for universal curves are hee = 0.100, h^o = 0.792, hoe = 
0.642, and Ko = 1.512. 

For Viola-Seaborg besides a = 1.3892, b = -13.862, c = 0.1086, d = 41.458 we introduce 
the additive hindrance constants CV of —0.073 for e-e, 0.484 for e-o, 0.478 for o-e, and 1.078 
for 0-0 nuclei. 

TABLE IL Standard deviations for semiempirical formula, and universal curves in the region of 
transuranium alpha emitters including superheavies {Z = 92 — 118). 



n parity 






<7 semFIS 


47 e-e 


0.282 


0.267 


0.164 


45 e-o 


0.562 


0.554 


0.540 


25 o-e 


0.608 


0.543 


0.514 


25 0-0 


0.597 


0.456 


0.492 



From the table |n] one can see how the universal curves reproduce better than the Viola- 
Seaborg formula the experimental data. A similar property is exhibited by the semFIS 
formula, which behaves very well for even-even alpha emitters. Moreover, the dependence 
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on the proton and neutron magic numbers of the semiempirical formula may be exploited to 
obtain informations about the values of the magic numbers which are not well known until 
now. 

In conclusion the new version of the SemFIS formula, having a smaller number of fitting 
parameters, gives slightly better results than the universal curve, both for a large number 
(493) of alpha emitters in the whole nuclear chart (Z = 52 — 118) and for 142 transuranium 
nuclei including superheavies {Z = 92 — 118). 
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FIGURE LEGEND 

Figure 1 Experimental half-lives for a-decay of 162 even-even (E-E), 122 even-odd 
(&0), 115 odd-even, and of 94 odd-odd (0-0) nuclides versus the neutron number of the 
daughter nucleus. The vertical bars correspond to spherical and deformed neutron magic 
numbers of the daughter nuclei A^^^ = 50, 82, 126, 152, 162, 172. They span a wide range of 
values between 10~^ and 10^^ s. The particularly strong shell effect at A'^^ = 126 is very 
clearly seen. 

Figure 2 Universal curves for a-decay of the same nuclides as in figure ^ in four groups 
of even-even, even-odd, odd-even and odd-odd nuclei. 

Figure 3 The deviations of a-decay half-lives calculated with the SemFIS semiem- 

pirical formula from the experimental values for even-even (top) and odd-even (bottom) 
nuclei. The vertical bars correspond to spherical and deformed neutron magic numbers of 
the daughter nuclei Nd = 50,82,126,152,162,172. There are nolonger large systematical 
discrepancies around magic numbers, because shell effects were taken into account by Sem- 
FIS. Calculations are performed with the new constants adjusted to fit the data of even-even 
nuclei. 
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